Role of tensor interaction in He isotopes with tensor-optimized shell model 
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We studied the role of the tensor interaction in He isotopes systematically on the basis of the 
tensor-optimized shell model (TOSM). We use a bare nucleon-nucleon interaction AV8' obtained 
from nucleon-nucleon scattering data. The short-range correlation is treated in the unitary corre- 
lation operator method (UCOM). Using the TOSM+UCOM approach, we investigate the role of 
tensor interaction on each spectrum in He isotopes. It is found that the tensor interaction enhances 
the LS splitting energy observed in ^He, in which the pi/2 and P3/2 orbits play different roles on the 
tensor correlation. In ^'^'*He, the low- lying states containing extra neutrons in the Pz/2 orbit gain 
the tensor contribution. On the other hand, the excited states containing extra neutrons in the pi/2 
orbit lose the tensor contribution due to the Pauli-blocking effect with the 2p2h states in the *He 
core configuration. 

PACS numbers: 21.60.Cs, 21.10.-k, 27.10. -|-h 27.20. H-n 
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INTRODUCTION 



It is an important subject in nuclear physics to un- 
derstand the nuclear structure from the viewpoint of the 
nucleon-nucleon (TV TV) interaction. The iViV interaction 
has distinctive features. There exist strong tensor inter- 
actions at long and intermediate distances caused by pion 
exchange and strong central repulsions at a short distance 
caused by quark dynamics. It is important to investigate 
the nuclear structure in relation to the above characteris- 
tics of the NN interaction [l|, [3] • Recently, it has become 
possible to calculate nuclei up to mass around A ^ \2 us- 
ing a interaction with the Green's Function Monte 
Carlo method (GFMC) [1, Q- K is, however, extremely 
time consuming to apply this method to heavier nuclei. 
It is desired to develop a new method to calculate nuclear 
structure with large nucleon numbers by taking care of 
the characteristic features of the A^A^ interaction. 

The presence of the tensor interaction in the NN in- 
teraction causes an explicit d-wave component in the rel- 
ative wave function in a nucleus, in particular, for the 
proton-neutron (pn) pair as seen in the dcuteron. The 
d-wave component in the deuteron is essential to bind 
the system via the sd coupling of the tensor interaction. 
This d-wavc component is found to be spatially compact 
as compared with the s-wave component due to the large 
momentum components brought by the tensor interac- 
tion 0]. This effect originates from the pseudo-scalar 
nature of the one-pion exchange. We mention also that 
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a large fraction of the pn pair is observed experimentally 
than pp or nn pairs in light nuclei 0] . This enhance- 
ment of the pn pair is hard to reproduce theoretically in a 
simple shell model (mean-field picture) M except for the 
rigorous method as the one of GFMC which treats 
the tensor interaction explicitly. It is important to make 
efforts to study the high momentum components caused 
by the tensor interaction in finite nuclei |9|. 

There are two important developments to proceed nu- 
clear structure calculations to heavy nuclei and to see the 
dynamics induced by the bare NN interaction. One is to 
find out that the strong tensor interaction is of interme- 
diate range and we are able to express the tensor correla- 
tion in a reasonable shell model space [13, [lli • We name 
this method as Tensor Optimized SheU Model (TOSM). 
The other is the Unitary Correlation Operator Method 
(UCOM) to treat the short-r ang e correlation caused by 
the short-range repulsion |12h14| . We shall combine two 
methods, TOSM and UCOM, to describe nuclei using 
bare A^A^ interaction. In the TOSM part, the wave func- 
tion is constructed in terms of the shell model basis states 
with full optimization of the two particle-two hole (2p2h) 
states. This means that there is no truncation of the par- 
ticle states within the 2p2h space of TOSM. In particular, 
the spatial shrinkage of the particle states is essential to 
obtain convergence of the tensor contribution involving 
the high momentum components [isl - fTsj . which is also 
seen in the deuteron. This treatment of the bare tensor 
interaction in TOSM corresponds to the one-pair approx- 
imation correlated by the tensor interaction jlSl Il9j . In a 
few-body viewpoint, Horii et al. have already confirmed 
the reliability of one-pair approximation to treat the bare 
tensor interaction using the few-body method for s-shell 
nuclei by reproducing more than 90% of the tensor cor- 
relation energy [2^. The explicit inclusion of the 2p2h 
states in the extended mean field model for heavy nuclei 
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has been formulated by Ogawa et al. [21|. 

So far, we have obtained successful results using TOSM 
for the investigation of the tensor correlations in He and 
Li isotopes. In ^He, we have shown the selectivity of 
(^1/2)^(51/2)^^ of 2p2h state with the pn pair induced 
by the tensor interaction and have recognized this corre- 
lation as deuteron-like one In ^He, this pn tensor 
correlation in ^He blocks the p -\ /■?, occupation of a last 
neutron by the Pauli-principle |22| . This blocking oc- 
curred in l/2~ state of ^He produces the p-wave splitting 
energy in ^He. In ^*^Li and ^^Li, we have performed the 
coupled two- and three-body analyses with the ^Li core 
described in TOSM, respectively. Similar to the case of 
^He, the pi/2 occupation of last one and two neutrons are 
blocked by the tensor correlation in ^Li. As a result, we 
have naturally explained the virtual s-state in lOLi, and 
the large s-wave mixing and a neutron halo formation in 
^^Li. We have also reproduced the various observables in 
"Li and "Li 

For the structures of He isotopes, many experiments 
have been reported in neutron-rich side |25l427| . How- 
ever, there are still contradictions in the observed energy 
levels and their spin assignments. We have performed 
the analyses of He isotopes on the basis of the cluster 
model treating the continuum states for extra neutrons 
j28l - l30j . We have reproduced various experimental data 
and predicted energy levels and decay widths. In these 
analyses, the ^He core is assumed to be the (Os)'^ configu- 
ration and the effects of 2p2h excitations are not consid- 
ered explicitly. In this approach, the dynamical coupling 
between the strong 2p2h excitations induced by the ten- 
sor interaction in ^He and the motions of extra neutrons 
in neutron-rich He isotopes cannot be treated. It is an 
interesting problem to take into account this coupling 
explicitly and to see how the coupling affects the energy 
levels and their tensor contributions in He isotopes. 

In this paper, we proceed with our study of the 
TOSM-I-UCOM approach to neutron-rich He isotopes 
and discuss their structures focusing on the roles of the 
tensor interaction on the energies and configurations. In 
the previous studies of '^'^He |22l. Isij. we have discussed 
the Pauli-blocking effect on the spectra of ^'^He based on 
the "^Hc+n+n three-body model, in which the ^He core is 
described in the TOSM. Semi-microscopic potentials be- 
tween ^He and an extra neutron was used, in which the 
tensor correlation is renormalized into the potential. In 
those analyses, the importance of the (pi/2)^ configura- 
tion of a pn pair in the ''He core is shown to produce the 
splitting energies between the states corresponding to the 
LS'-partners in ^'^He. In the present study, we perform 
a full microscopic analysis by using the TOSM+UCOM 
approach with the bare NN interaction for He isotopes 
and see how the present method works. We discuss the 
role of the pi/2 orbit excited by the tensor interaction 
and its structure dependence in the individual states ob- 
tained in ^^®He. We also see the kinetic energy behavior 
and the high momentum component in each state in re- 
lation with the tensor correlation. This subject becomes 



a foundation of the nuclear structure study on the basis 
of TOSM+UCOM. 

In Sec. ini we explain the method of the 
TOSM-hUCOM approach for He isotopes. In SecHIl we 
show the level structures of He isotopes and discuss their 
characteristics in relation with the tensor interaction. A 
summary is given in Sec. IIVI 

II. MODEL 

A. Tensor-optimized shell model (TOSM) for He 
isotopes 

We explain the framework of TOSM. We shall begin 
with a many-body Hamiltonian, 

H = 5^T,-re,„.+^F,„ (1) 

v., = vf^+v^+v^ + v^^-K (2) 

Here, T,; is the kinetic energy of each nucleon with Tc.m. 

being the center of mass kinetic energy. We take a bare 

interaction Vij such as AV8' [3| consisting of central 

tensor vj, and spin-orbit vj"f terms and w,^'™'' the 
''J V 

Coulomb term. We obtain the many-body wave function 
with the Schrodinger equation = E'i/. 
We give the TOSM wave function \1/ as 

+ Y,Akj2p2h;k2), (3) 

\lplh;k,) = |Cj(^i)xC?(ri)), (4) 
\2p2h;k2) = |V'Sj(n)VS^(r2)xC3^(ri)Vi^:(r2)).(5) 

Here, we omit the symbols of angular momentum cou- 
pling and total isospin for simplicity. The states 
OpO/i; fco) are the OpOh shell model states. The config- 
urations \lplh; ki) and \2p2h; ^2) are the Iplh and 2p2h 
states with various radial components for particle states. 
The vectors ri and r2 represent the positions of nucleons. 
The labels fcoj ^i and fc2 are the representative indices to 
distinguish various configurations. We take the available 
configurations with the fixed spin and parity J'" and the 
isospin of the total wave function ^P. The amplitudes 
{AfcQ, Afcj , Afcg} of the configurations are variational co- 
efficients. The basis functions ip2 and V'" are to describe 
particle and hole states of one nucleon with quantum 
numbers n and a, respectively; the index n is to distin- 
guish different radial basis functions and a are the sets of 
the spin-isospin quantum numbers to distinguish the or- 
bits. In TOSM, the hole states are described by harmonic 
oscillator basis states, whose length parameters are de- 
termined independently for each orbit in each J"^ state so 
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as to minimize the total energy of the system. For ^He, 
the OpOh configuration is given by the (Os)^ wave func- 
tion and for ^~*He, the OpOh configurations are given in 
the Os + Op space with Ohuj excitations. The available 
configurations for particle states are taken into account 
within the Iplh and 2p2h states with up to the partial 
wave Lmax of VSj which determines the model space of 
TOSM. We check the convergence of the solutions by in- 
creasing i,nax- The 2p2h states play an important role on 
the description of the strong tensor correlation in TOSM 
and the Iplh states can improve the radial properties of 
the hole states. 

Wc employ the Gaussian expansion method to de- 
scribe single-particle basis states tp^ particle states 
in TOSM [33 . [3^. Each Gaussian basis function has 
the form of a nodeless harmonic oscillator wave functions 
(HOWF), except for the Is and Ip orbits. When we su- 
perpose a sufficient number of Gaussian basis functions 
with various length parameters, the radial component of 
the particle states can be fully optimized in every con- 
figuration of TOSM with respect to the Hamiltonian in 
Eq. ([T]). The basis states for particle states should be or- 
thogonal to the hole states having HOWF. This condition 
is imposed by using the Gram-Schmidt orthonormaliza- 
tion which is explained later. Technically, in order to use 
the non-orthogonal Gaussian basis function in the shell 
model framework, we construct the following orthonor- 
malized single-particle basis function used in Eqs. ^ 
and ([5|) by using a linear combination of Gaussian bases 
{(f>a} with length parameter 6q,.^. 



(CIV'S') = Sn,n'Sa,a', 

for n = 1 , • ■ • , Na , 



(6) 
(7) 



where Na is a number of basis functions for the orbit 
a, and v is an index that distinguishes the bases with 
different values olba^^. The explicit form of the Gaussian 
basis function is expressed as 



2 h-) 



(2i+3) 



r(/ + 3/2) 



(9) 



where I and j are the orbital and total angular momenta 
of the basis states, respectively, and tz is the projection 
of the nucleon isospin. The weight coefficients {d" j,} are 
determined to satisfy the overlap condition in Eq. 
This is done by solving the eigenvalue problem of the 
norm matrix of the Gaussian basis set in Eq. ^ with 
the dimension Na- Another method is the Gram-Schmidt 
orthonormalization, in which the basis function -02 is ex- 
panded by the Gaussian basis functions with the number 
n These two methods give the different {d^ ^} to 

prepare {V'Sli but, are equivalent to obtaining the vari- 
ational solutions of the total wave function vj/ in Eq. ^ . 



This property is satisfied for any Gaussian basis number 
Na- Following this procedure, we obtain the new single- 
particle basis states {i/'Sl in Eq. ^ used in TOSM. The 
particle states in each configuration are determined inde- 
pendently in TOSM from the variational principle. 

Wc should construct Gaussian basis functions of par- 
ticle states to be orthogonal to the occupied states, the 
Os and the Op orbits in He isotopes. We explain this 
procedure as follows; For the Os orbit in the hole state, 
we employ one Gaussian basis function, namely, HOWF 
with length 60s, 1^=1 = b^s- For the lsj^/2 basis states in 
particle states, we introduce an extended Is basis func- 
tion orthogonal to the Osi/2 state and that possesses a 
length parameter his^u that can differ from 6os- In the ex- 
tended Is basis functions, the polynomial part is changed 
from the usual Is basis states to satisfy the conditions of 
the normalization and the orthogonality to the Os state 
p^ . For the Ip states, wc take the same method as 
used for the Is case. One-body and two-body matrix 
elements in the Hamiltonian arc analytically calculated 
using the Gaussian bases, whose explicit forms are given 
in appendix of Ref. [Tl| for central, LS and tensor in- 
teractions, respectively. In the numerical calculation, we 
prepare 10 Gaussian basis functions at most with various 
range parameters to get a convergence of the energy and 
Hamiltonian components. 

Wc furthermore take care of the center-of-mass excita- 
tions. Toward this end, we take the well-tested method 
of introducing a Hamiltonian of center-of-mass motion 
in the many-body Hamiltonian known as the Lawson 
method [13]. In the present study, we take the value 
of fiLu for the center of mass motion as the averaged one 
used in the Os and Op orbits in the OpO/i states with the 
weight of the occupation numbers in each orbit. We have 
checked that this choice of Tiuj gives variationally better 
solutions in the present scheme. Adding this center of 
mass Hamiltonian as the Lagrange multiplier to the orig- 
inal Hamiltonian in Eq. ([T]) , we can effectively project out 
only the lowest HO state for the center-of-mass motion. 
The value of the Lagrange multiplier is taken typically 
as 10-20 to suppress the center of mass excitations below 
100 keV. 

The variation of the energy expectation value with re- 
spect to the total wave function in Eq. ([3]) is given 
by 



which leads to the following equations: 



ki 







= for i = 0,1,2. 



(10) 

(11) 
(12) 



Here, total energy is represented by E. The parameters 
{6q,i/} for the Gaussian bases appear in non-linear forms 
in the energy expectation value. We solve two kinds of 
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variational equations in Eqs. pT|) and ([T^ in the follow- 
ing steps. First, fixing all the length parameters ba^u and 
the partial waves of basis states up to -Lmax, we solve the 
linear equation for {Ak- } as an eigenvalue problem for H. 
We thereby obtain the eigenvalue E, which is functions 
of {ba.iy} and £max- Ncxt, we try to adopt various sets of 
the length parameters {ba,u} and expand the configura- 
tion space of TOSM by increasing imax in order to find 
the solution which minimizes the total energy E. In the 
TOSM wave function, wc can describe the spatial shrink- 
age of particle states with an appropriate radial form in 
each configuration within 2p2h space, which is important 
to describe the tensor correlation [10|. 



B. Unitary Correlation Operator Method (UCOM) 

We explain UCOM for the short-range central correla- 
tion |13 - [l4| . in which the following unitary operator C 
is introduced 



C = exp(-i^g„) 



(13) 



We express the correlated wave function in terms of 
less sophisticated wave function $ as ^ = C$. Hence, 
the Schrodinger equation becomes = E^ where the 
transformed Hamiltonian is given as iJ = C'^HC. Since 
the operator C is expressed with a two-body operator 
in the exponential, it is in principle a many-body opera- 
tor. In case of the short-range correlation, we are able to 
truncate the modified operators at the level of two-body 
operators [l^ . 

Two-body Hermite operator g in Eq. (|13p is defined as 



9 = \ {Prs{r) + s{r)pr] 



dR+{r) s{R+{r)) 



dr 



s{r) 



(14) 



(15) 



where the operator pr is the radial component of the 
relative momentum and is conjugate to the relative co- 
ordinate r. The function s{r) expresses the amount of 
the shift of the relative wave function which depends on 
the relative coordinate. In UCOM, we use R+{r) in- 
stead of s(r), where the function R+{r) corresponds to 
the transformed relative coordinate from the original r. 
This i?+ (r) represents the correlation function to reduce 
the short-range amplitude of the relative wave functions 
so as to avoid the short-range repulsion in the NN in- 
teraction. The explicit transformation of the operator is 
given in Refs. [H, [l^. We use the TOSM basis states to 
describe $ which includes the tensor correlation. 

In UCOM, the shift operator g in Eq. is intro- 
duced for every nuclcon pair in nuclei. The amount of 
the shifts, namely the the function R+{r) is determined 
variationally. We parametrize i?+ (r) in the same manner 
as proposed by Feldmeier and Neff [13, fisj . We assume 



TABLE I: Optimized parameters in R+{r) 
TOSM+UCOM in four spin-isospin channels. 



used for 





a 


13 


7 


singlet even 


1.32 


0.88 


0.36 


triplet even 


1.33 


0.93 


0.41 


singlet odd 


1.57 


1.26 


0.73 


triplet odd 


1.18 


1.39 


0.53 




2.5 



FIG. 1: The UCOM functions R+{r) used in even C'''^E) and 
odd (^'■'O) channels. The differences between R+(r) and the 
original coordinate r are displayed. 



the following forms for even and odd channels, respec- 
tively. 

i?-on(^) ^ ^ + ^(^^ycxphexp(^)], (16) 
R°+'^{r) = r-Hafl-exp(--))exp[-exp(^)](17) 
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Here, a, /3, 7 are the variational parameters to optimize 
the R+{r) functions and minimize the total energy of 
the system. They are independently determined for four 
channels of spin-isospin pair. In the actual procedure of 
the variation, once we fix the parameters in R^{r), we 
solve the eigenvalue problem of the transformed Hamilto- 
nian using Eqs. ([TT]) and ([T^ and determine the config- 
uration mixing in TOSM. Next, we try to search various 
sets of the R+{r) parameters to minimize the total en- 
ergy. 

In the framework of UCOM, it is generally possible 
to introduce the partial wave dependence in R+{r) and 
then the i?+ (r) functions are determined in each relative 
partial wave of the nucleon pair. In the previous study 
UJ, we have performed the extension of UCOM by tak- 
ing care of the characteristics of the short range correla- 
tion. One of the simplest extensions is UCOM for only 
s-wave relative motion. The partial waves with finite or- 
bital angular momentum / except for the s-wave {I = 0) 
have r' behavior of the relative wave function near the 
origin. This r' behavior cuts down the effect of the short- 
range hard core in the NN interaction largely. On the 
other hand, only the s-wave function is finite at the ori- 
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gin and its behavior near the origin is determined by the 
hard core dynamics. Considering this physical charac- 
teristics of the s-wave, we have introduced 5-UCOM, in 
which the UCOM transformation is performed for only 
the relative s-wave component in the even channel. In 
the ordinary UCOM, there is too large a removal of the 
short-range part of the relative wave functions, in par- 
ticular, in the c?-wave part used in the sd coupling of 
the tensor matrix elements, where the tensor interaction 
possesses some amount of the strength. In order to over- 
come this feature, we have newly introduced S'-UCOM 
and confirmed that the situation was improved for the 
binding energy by a few MeV in "He [ll|. In TOSM 
with S'-UCOM, we can optimize the R+{r) functions and 
the optimized three parameters in S'-UCOM for ^He are 
listed in Table IT] The demonstration of the calculation to 
search for the energy minimum is shown in Ref. . The 
adopted R+{r) in S-UCOM for four channels are shown 
in Fig. [T] In the present analysis, we use these R+{r) 
functions for every states of He isotopes. To simplify the 
numerical calculation, we adopt the ordinary UCOM for 
the central correlation part instead of the S-UCOM in 
this analysis. 



III. RESULTS 
A. *He 

We show first the resuhs of ^He using TOSM-hUCOM. 
We then explain the structures of '''~®He and discuss how 
their structures change in comparison with that of ^He. 
We explain numerical results of ^He using the AV8' in- 
teraction, which consists of central, LS and tensor terms 
and is used in the calculation given by Kamada et al, 
where the Coulomb term is ignored 0. In Fig. 1 the 
convergence of the energy of ''^He is shown as a function of 
^max, the maximum orbital angular momentum of par- 
ticle states in TOSM. We get a good convergence for 
^Hc. The Hamiltonian components are shown in Table 
im in comparison with the stochastic variational method 
(SVM) using correlated Gaussian basis functions [ssl. [36j . 
which is one of the rigorous calculations ■ The matter 
radius of "^He is obtained as 1.52 fm in TOSM+UCOM. 
When we apply the S-UCOM instead of ordinary UCOM, 
the energy gain is about 2 MeV in total energy and gets 
closer to the rigorous value [ll|. In particular, contri- 
bution from the tensor interaction becomes large due to 
the improvement of the sd coupling of the tensor matrix 
elements as was mentioned. 

As for the comparison with the rigorous calculation, 
we see that contribution from the central interaction 
in TOSM+UCOM/S-UCOM satisfies the rigorous value. 
On the other hand, the tensor energy and the kinetic en- 
ergy show some shortage from the rigorous values even 
in the S-UCOM case. One of the possible explainations 
for the shortage of those energies in this approach is the 
treatment of the short-range part of the tensor correla- 



TABLE II: Various energy components in MeV for *He. 





Energy 


Kinetic Central 


Tensor 


LS 


TOSM-fUCOM 
TOSM-hS-UCOM 
SVM \2} 


-20.46 
-22.30 
-25.92 


86.95 -54.63 
90.50 -55.71 
102.35 -55.23 


-51.06 
-54.55 
-68.32 


-1.73 
-2.53 
-4.71 



tions. Although the dominant part of the tensor inter- 
action is of intermediate and long ranges, there remains 
small strength in the short-range part of the tensor in- 
teraction, which can couple with the short-range correla- 
tions. Horii et al. estimate the amount of this coupling 
by using the few-body method of SVM [g^l, in which 
the short-range repulsion and the tensor interaction are 
directly treated. In Ref. [1^, they propose the one- 
pair approximation of tensor coupling in the few-body 
wave function, named tensor-optimized few-body model 
(TOFM). It is shown that the TOFM gives the good 
binding energy of *He of 24.05 MeV with AV8' as com- 
pared with the rigorous calculation. The physical con- 
cept of TOFM is the same as that of TOSM except for 
the use of UCOM. These results imply that the UCOM 
transformation for short-range correlation makes the en- 
ergy loss of about 1.5 MeV in '''He in comparison with 
the TOFM result, in particular, for the tensor contri- 
bution. The remaining shortage of the binding energy 
with respect to the rigorous calculation shown in Table 
nil should come from higher configurations beyond 2p2h 
excitations in TOSM. Horii et al. also confirmed that the 
two-pair tensor coupling in *He by extending TOFM, re- 
produces the rigorous results within 200 keV [s^, which 
corresponds to the 4p4/i mixing in TOSM. 

Considering the difference between the results of 
TOSM and TOFM, three-body term of the UCOM trans- 
formation is one of the possibilities to overcome the lack 
of energy from UCOM [ll|, [T^. This three-body term 
may contribute to the increase the tensor energy and 
bring more high momentum components in the wave 
function. It is also noted that it has been discussed 
that the three-body term of UCOM can work repulsively 
more or less for the Hamiltonian having only the central 
NN interaction [l^] and that a similar trend can be seen 
in UCOM including the tensor-type transformation [l^ . 
Another possibility is the improvement of the correlation 
function R^(r). In this calculation, the functional forms 
of R+{r) shown in Eqs. and p7)) are introduced 

from the consideration of the short-range behavior of the 
two-body system only with the central TV TV interaction 
[T^ . It would be interesting to study an appropriate form 
of i?-)- (r) which is suitable for the Hamiltonian including 
the tensor interaction explicitly, and discuss the effect of 
the many-body term of UCOM. The detailed analysis of 
the result obtained in TOFM [g^l gives some ideas to de- 
termine the form of R+{r) from the variational principle. 

To see the properties of the '*He wave function, dom- 
inant configurations are listed in Table IIIII with their 
probabilities. For the particle states, the various radial 
components are summed up for each orbit. It is found 
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TABLE IV: Occupation numbers in each orbit in *He. 



10 12 14 



FIG. 2: Energy of *He in TOSM+UCOM as a function of the 
maximum orbital angular momentum Lmax- 



TABLE III: Mixing probabilities in % for "^He with 
TOSM+UCOM, where the two subscripts 00 and 10 are the 
spin-isospin quantum numbers. 



(0^0 
(0s)w(0Pi/2)?o 

(Os)-o^[(lSl/2)(Od3/2)]lO 
(Os)ro'[(Op3/2)(0/5/2)]lo 
(0s)r<f[(0Pl/2)(0p3/2)ll0 
(Os)-(f[(0d5/2)(0g7/2)]l0 



remaining part 



84.14 
2.32 
2.20 
1.82 
1.21 
0.78 



7.53 



that the specific 2p2h states such as (0s)]]3^(0pi/2)io show 
large probabilities and these configurations are essential 
to produce the tensor correlation in ^He because of the 
coupling by the tensor operator [lol [l6j . These dominant 
2p2h states commonly have the isoscalar nature and cor- 
respond to the excitations of a pn pair, which indicates 
the deuteron-like correlation. It has already been shown 
that the particle states in those 2p2h states have the spa- 
tially compact radial wave functions [l^jIUIl^- I* is also 
found that the QpQh-2p2h coupling exhausts 94% of the 
tensor energy. Occupation numbers in various nuclcon 
orbits are shown in Table IIVI The pi/2 orbit has the 
largest occupation number among particle states because 
of the large 2p2h mixings including the component. 
This feature of 2p2h excitations plays an important role 
to determine the structures of ^^^He as will be discussed 
later. 



B. Energy spectra of He isotopes 

We explain the results of neutron-rich He isotopes us- 
ing TOSM-fUCOM with the AV8' interaction, where 
Lmax is taken as 10 to get a sufficient convergence. We 
show the energy spectra with respect to the ^Hc ground 
state energy in Fig. [3] The excitation energies of ^~^Hc 
are shown in Fig. 21 The matter radii of ^He and ®He are 
listed in Table |Vl which are slightly smaller than those 
of the experiments [ssl - lioj . We also show the results of 



*He(J") 


Sl/2 


Pl/2 


P3/2 


^3/2 


^5/2 


0+ 


3.80 


0.06 


0.04 


0.04 


0.01 



TABLE V: Matter radii of '^He and '^He in comparison with 
the cluster orbital shell model (COSM) [s^] and the experi- 
ments; afSal, b[391, cliyl. Units are in fm. 





Present 


COSM 


Experiment 




«He 


2.27 


2.37 


2.33(4)=^ 2.45(10)'' 


2.37(5)=^ 


^He 


2.44 


2.52 


2.49(4)=^ 2.53(8)'' 


2.49(4)"= 



cluster orbital shell model (COSM) [33|, in which the 
spatial extension of extra neutrons are fully described. 
This point will be discussed later. For charge radius, 
the precise data of ^'^He arc reported [4l|. To obtain 
the charge radius it is necessary to calculate the isospin- 
breaking matrix elements, which are not yet obtained in 
the present TOSM. 

It is found that our results underestimate the binding 
energies of He isotopes and its amount becomes larger for 
neutron-rich side. On the other hand, for excitation en- 
ergies shown in Fig. 21 we sec a good correspondence to 
the experimental energy levels. This result indicates that 
the level spacing is described well in TOSM-I-UCOM, so 
that we can discuss the structure differences between en- 
ergy levels. Before doing this, we discuss several rea- 
sons of under-binding, which can contribute to the bulk 
part of the binding energies, such as the continuum effect 
of extra neutrons having a weak binding nature, higher 
configurations beyond 2p2h states in TOSM, the genuine 
three-body interaction, and the improvement of UCOM. 

As for the continuum effect, the spatial extension of 
extra neutrons is explicitly treated in the TOSM ba- 
sis within the 2p2h excitations using the Gaussian ex- 
pansion method by taking long-range parameters in Eq. 

However, the continuum effect of extra neutrons is 
not sufficiently described in the TOSM approximation. 
When the 2p2h excitations in TOSM are used to describe 
the pn tensor correlation in ^He, the extra neutrons still 
occupy the p-orbits in the hole states. This means the 
continuum effect is not taken into account together with 
the large energy gain due to the tensor correlation in 
TOSM. Hence, we estimate the correlation energy from 
the continuum effect by using COSM [s^], in which each 
motion of extra neutrons is solved exactly around the 
core with no tensor correlation, and the continuum effect 
of multi-neutrons is fully included. It is obtained that 
the correlation energies beyond the p-shell configuration 
of extra neutrons are estimated as 3.3 MeV for ^He and 
6.6 MeV for ^He in their ground states, respectively [30| . 

It is important then to estimate how many particle- 
many hole configurations arc necessary to include the 
continuum effect of extra neutrons. For ^He, when we 
limit COSM up to 2p2h excitations of extra four neu- 
trons, the energy gain is 5.9 MeV, which is about 90% 
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of the correlation energy of the continuum effect of 6.6 
MeV. This result indicates that a large fraction of the 
continuum effect can be described within the 2p2h space 
from the p-shell neutrons. This fact suggests that the 
extension of TOSM to include 4p4/i configurations as a 
whole can explain this feature of extra neutrons for He 
isotopes with two-kinds of 2p2h contributions; a pn pair 
with r = in the '^He core and a nn pair with T = 1 in 
the p-shcU neutrons. Furthermore, to satisfy the bound- 
ary condition of neutron emissions, it is necessary to 
connect the TOSM basis to the correct asymptotic wave 
functions such as the ^He-f n-t-n three-body state for ^He 
[sot . This is important to describe the tail behavior of 
the neutron wave functions. 

From the GFMC calculation we can simply esti- 
mate the effect of the genuine three-body interaction on 
the binding energy which increases with the mass num- 
ber. For *He, for example, the relative energy gain from 
three-body interaction is about 5.5 MeV as compared 
with ''He. This result indicates that for ^He, the total 
energy gain of about 12 MeV is expected considering two 
effects as the continuum states and the genuine three- 
body interaction. Hence, the remaining missing energy 
from the experimental value is estimated to be as small 
as 3-4 MeV for **Hc from Fig. [3l 

For the UCOM part, wc arc able to increase the ten- 
sor contributions and the binding energies by using S- 
UCOM, which can improve short-range part of the sd 
coupling of the tensor interaction. The improvement of 
the UCOM correlation function R+{r) provides another 
possible way to increase the energy, according to the re- 
sults obtained in TOFM as explained above. 

In He isotopes, most of the states are observed as res- 
onances and in this analysis we describe those states 
within the bound state approximation. This treatment 
makes it possible to focus the discussion on the inter- 
nal structures of each level in He isotopes systematically. 
In the previous studies [H, [H, [U , we include the con- 
tinuum effects explicitly, using the TOSM wave func- 
tion as the core nucleus in the extended two- and three- 
body cluster models. In those analyses, we have success- 
fully explained the ^Hc+?i scattering phase shifts, neu- 
tron halo formation and three-body Coulomb breakup 
strengths in halo nuclei ®He and ^^Li and so on, al- 
though the systems are not described starting from the 
single NN interaction like the present TOSM-I-UCOM. 
It was found that the 2p2h components induced by the 
tensor interaction play an important role in reproducing 
the physical observables in those studies. Furthermore, 
the resulting characteristics of the tensor interaction are 
essentially the same as that obtained in the present anal- 
ysis. Based on these achievements, in the present study 
it is important to understand the role of tensor interac- 
tion in each state of He isotopes systematically fully using 
TOSM-f UCOM starting from the bare NN interaction. 
The present approach can also be extended to the anal- 
ysis which includes the continuum states explicitly, such 
as the scattering state in the ''He-l-n system [i^ ]. 
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FIG. 3: Energy levels of He isotopes using AV8', normalized 
to the ^He ground state energy. Experimental data are taken 
from Refs. BlSl- 
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FIG. 5: Energy levels of the He isotopes using the Minnesota 
interaction normalized to the *He ground state energy. 



To see the interaction dependence, the energy spectra 
using the Minnesota (MN) interaction are shown in Fig. 
[5] The MN two-body interaction consists of the cen- 
tral and LS parts with u parameter being 0.95 [43l |4^. 
We use Set III of the MN interaction corresponding to 
u = 0.95 for the LS forceQ. There is no explicit tensor 
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correlation in the MN interaction. We do not use UCOM 
for the calculation using the MN interaction, because the 
MN interaction does not have a strong short-range re- 
pulsion. It was already confirmed that TOSM can repro- 
duce the binding energy and the radius of "^Hc obtained 
by few-body SVM approach with a central interaction 
case [3^. The difference of binding energy with MN is 
about 200 keV in comparison with the SVM approach 
[ssj . Using the MN interaction, the binding energy of 
^He is obtained as 30.55 McV without the Coulomb in- 
teraction. The radius of *He is obtained as 1.39 fm, which 
is smaller than the experiment by about 0.1 fm. For ^He 
and '^He, their matter radii arc obtained as 1.93 fm and 
1.91 f m, respe ctively, which are smaller than the experi- 
ments [38l - l40l | shown in Table fVl bv about 0.5 fm. From 
the energy spectra shown in Fig. [Sj it is found that the 
reproduction of ^He energy is good, on the other hand 
the energies of ^He and ®He are overestimated. One of 
the reasons is the too large splitting energy between P3/2 
and Pi/2 components mainly given by the LS interaction. 
When wc adjust the MN interaction such as changing the 
u parameter and the strength of the LS interaction, it 
is still difficult to reproduce the whole trend of the en- 
ergies of the ground and excited states of He isotopes 
consistently. One of the reasons to explain this tendency 
is that the MN interaction is originally formulated as an 
effective interaction when "^He is assumed to be the (Os)'* 
configuration [i^ l45l - l47j . This condition is not imposed 
in TOSM. Another difficulty of the MN interaction is 
the radius property, which is related to the nuclear sat- 
uration density. The MN interaction is known to give 
spatially compact wave functions than the experimental 
values, even for the (O.s)'* state of '^He (45j . 



C. Hamiltonian components 

We discuss the structures of each level of He isotopes. 
Various energy components in ^"^He are listed in Tables 
ED Eni rvini and HXl which are useful to discuss the 
structures of individual states. In this calculation, we 
take the common length parameters of the hole states, 
Os and Qp orbits. This is done to exclude the continuum 
effect, which produces a few MeV energy gain in ^He as 
shown in Fig[31 and to focus our discussion on the internal 
structures of He isotopes. We take 1.5 fm of the length 
parameter of the hole states, which is the optimized value 
of the Os orbit in ^He. This length parameter corresponds 
to huj as 18.43 MeV. In the full calculations explained 
above, we minimize the energies by optimizing the length 
parameters of all the hole states for each J'^ state of He 
isotopes. 

In Table I VII we compare various energy components 
in 3/2^ and 1/2" states of ^He measured from those of 
"'He. The LS splitting energy obtained is about 3 MeV. 
We discuss the effect of the tensor interaction on this 
splitting energy between the two states. A large differ- 
ence is seen in the tensor energy of the 3/2 ^ state as 



TABLE VI: Various energy components in ''He measured from 
those of the ''He ground state. Energy units are given in MeV 
and fiu — 18.43 MeV. Dominant configurations of the extra 
neutron are shown also. 



5He(J") 


Config. 


Energy Kinetic Central Tensor 


LS 


3/2- 
1/2- 


P3/2 

Pl/2 


6.97 24.14 -8.99 -5.60 
10.05 17.53 -6.96 -1.11 


-2.58 
1.04 



TABLE VII: Various energy components in ^He measured 
from those of the ^He ground state. Units are in MeV. 



''He(J'') 


Config. 


Energy 


Kinetic Central 


Tensor 


LS 


0^ 


(P3/2)^ 


8.95 


53.04 


-27.75 


-12.02 


-4.04 


Ot 


(P1/2)' 


21.90 


34.30 


-14.06 


-0.17 


2.11 


1+ 


(P3/2)(Pl/2) 


17.29 


42.90 


-15.98 


-8.49 


-0.86 




(P3/2)^ 


11.40 


52.41 


-22.93 


-12.80 


-4.99 




(P3/2)(Pl/2) 


16.07 


45.06 


-18.06 


-8.54 


-2.10 



compared with that of 1/2^. The larger contribution of 
the tensor interaction in 3/2" brings the enhancement 
of the kinetic energy because of the involvement of high 
momentum components from the tensor interaction. The 
amount of the enhanced kinetic energy is 24 MeV, which 
is larger than huj. For 1/2-, on the other hand, the en- 
ergy gain from the tensor interaction is small and the 
enhancement of the kinetic energy is 17.5 MeV, which 
is close to the value of Tiuj. These results are related to 
the larger amoimt of the Pi/2 component in *He than 
the P3/2 one as shown in Table Hill This situation is ex- 
plained in Fig. [6l where the last neutron in ^He occupies 
the orbit in the left panel, this neutron does not dis- 
turb the ''He structure. Hence, the P3/2 occupied state 
gains an additional tensor energy without disturbing the 
large energy gain in '*Hc. On the other hand, in case of 
the dominant Pi/2 occupation for the 1/2- state shown 
in the right panel of Fig. [6l this neutron blocks some 
component of the spatially compact Pi/2 neutron in the 
2p2h excitations of the ''He core configuration because of 
the small degeneracy of the pi/2 orbit. This effect dy- 
namically produces the Pauli-blocking and reduces the 
total binding energy of ^He. As a result, the last neu- 
tron located in the pi/2 orbit should be orthogonal to the 
excited P\/2 orbit in *He and the tensor interaction does 
not gain the energy in ^He(l/2-). The coupling behavior 
between *He and a last neutron explains the difference in 
the Hamiltonian components in two states of ^He, which 
results in the LS splitting energy as a net value. 

In order to show directly the relation between the LS 
splitting energy in ^He and the tensor interaction, we 
show the LS splitting energy of ^He between 3/2^ and 
1/2- states in Fig. [7] by changing the strength of ten- 
sor interaction. In this process, we always optimize the 
length parameters of hole states independently in the 
same manner as is done in Fig. [31 When the tensor 
interaction is enhanced, the LS splitting energy becomes 
larger, because of the different couplings between ^He 
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FIG. 6: Configurations in ^He. Open circles with solid lines 
indicates the extra neutrons. 
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FIG. 7: LS splitting energy in ^He as a function of the 
strength of the tensor interaction Xtcnsor , which is multiplied 
to the tensor interaction in AV8'. 



and a last neutron occupied in the pi/2 or orbit. This 
result shows that the tensor interaction plays a decisive 
role to create the LS splitting energy in ^He [Zsl - fsoj . In 
our previous studies |22l. |3]|. we have estimated the con- 
tribution of this tensor correlation to the observed phase 
shifts of the ''He-n scattering, which produces about 30% 
of the splitting. A similar tendency is confirmed for other 
He isotopes, such as the level spacing between the and 
O2 states in ^He. 

For ^He, various energy components are shown in Ta- 
ble I VIII similarly to ^He. We can classify the obtained 
five states of ^He into three groups by seeing the ten- 
sor energies. The first group is the ground 0^ and 2^ 
states, where the extra two neutrons dominantly occupy 
the (p3/2)^ configuration. In this case the tensor corre- 
lation in ■*IIe is hardly disturbed due to a small mixing 
of the neutron in ^He, thus the tensor interaction 
acts attractively to the binding energy as illustrated in 
the left panel of Fig. [51 The difference of the central 
interaction energies in two states indicates the existence 
of the neutron pairing correlation in the 0^ state. In 
these states, we see also the enhancement of the kinetic 
energies coming from the high momentum component as- 
sociated with the tensor interaction beyond the value of 
2 X hw. The second group is the 1"*" and 2^ states, where 
the extra two neutrons occupy the (pi/2 P3/2) configu- 
ration. In these states, similarly to the ^IIe(l/2-) case. 
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FIG. 8: Configurations in He. Open circles with solid lines 
indicates the extra neutrons. 



TABLE VIII: Various energy components in '^He measured 
from those of the *He ground state. Units are in MeV. 

'He(J'^) Config. Energy Kinetic Central Tensor LS 



3/2r 


(P3/2)' 


14.97 


82.89 


-39.13 


-19.89 


-8.63 


3/22- 


(P3/2)i(Pl/2) 


21.44 


74.14 


-33.64 


-15.98 


-2.80 


3/23- 


(P3/2)'(Pl/2)' 


27.50 


65.82 


-27.09 


-9.96 


-0.98 


1/2- 


(P3/2)o(Pl/2) 


19.07 


76.81 


-38.52 


-15.74 


-3.20 


5/2^ 


(P3/2)2(Pl/2) 


19.92 


77.53 


-36.31 


-16.60 


-4.37 


TABLE IX: Various ener 


gy components in *He measured from 


those of the *He ground state. 


Units are 


in MeV. 




»He(J") 


Config. 


Energy Kinetic Central 


Tensor 


LS 


0^ 


(P3/2)^ 


17.57 


114.87 


-59.26 


-25.97 


-11.79 




(P3/2)^(Pl/2)^ 


29.50 


101.03 


-51.60 


-17.10 


-2.51 


1+ 


(P3/2)^(pi/2) 


25.08 


106.52 


-50.29 


-23.75 


-7.11 


2^ 


(P3/2)^(Pl/2) 


22.54 


109.60 


-56.22 


-23.42 


-7.13 


2J 


(P3/2)^bl/2)^ 


30.13 


100.15 


-49.27 


-18.08 


-2.38 



tensor energy and kinetic energy decrease from those in 
the first group because one of the extra neutrons occu- 
pies the Pi/2 orbit and the blocking of this orbit occurs. 
The third group is the state with the main occupa- 
tion probability in the configuration for the extra 
two neutrons. In this state, the Pauli-blocking effect to 
block the tensor correlation strongly occurs as shown in 
the right panel of Fig. [H As a result, the tensor interac- 
tion cannot produce the energy gain from the *Hc case 
and the kinetic energy shows the lowest value among the 
^He states. This blocking dynamically produces the LS 
splitting energy and explains the level spacing in ^He. 
From these results of ^He, it is found that the configura- 
tions of extra two neutrons in ^He are essential to explain 
the difference seen in various Hamiltonian components of 
the individual state shown in Table IVIll Enhancement of 
the kinetic energy as the increase of the high momentum 
component beyond the Tioj effect is also originated from 
the enhancement of the tensor energy. 

The similar effect about the roles of tensor correlation 
as seen in ^He and ®Hc is obtained for ^He and ®He. For 
^He, there exist three kinds of groups. The first group, 
the ground 3/2 state mostly exhausts the tensor en- 
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ergy due to the dominant (^3/2)"^ configuration of the 
extra three neutrons. The kinetic energy also shows the 
large value beyond 3 x hoj. The second group, the 3/2^, 
l/2~ and 5/2~ states show the similar energy compo- 
nents because they have commonly the main configura- 
tion of (^3/2)^(^1/2) of the extra three neutrons. The 
third group, the 3/2'^ state shows the smallest tensor 
and kinetic energies among the ^He states. This state 
mostly has the Pi/2 component in ^He. For ®He, five 
states are obtained and the three kinds of groups on the 
contributions of tensor interaction are seen. The ground 
0^ state has mainly the (^3/2)^ configuration of extra 
four neutrons. Due to this configuration, the 0^ state 
mostly exhausts the tensor energy and simultaneously a 
large kinetic energy beyond that of 4 x hcj. The 1"*" and 
2^ states form the second group having mainly three P3/2 
orbit neutrons. The Oj and 2 J states show small tensor 
and kinetic energies due to the large Pi/2 component. 

As conclusions of the roles of the tensor interaction in 
He isotopes, the deuteron-like pn pair of the tensor corre- 
lation in ^He is coupled with the motion of extra neutrons 
in neutron-rich He isotopes. There are two-kinds of cou- 
plings. One is when the extra neutrons occupy the P3/2 
orbit, these neutrons do not disturb the tensor correlation 
in ^He due to the small occupation of the P3/2 neutron 
orbit in '*He. As a result, the total system can enhance 
the tensor contribution and simultaneously the high mo- 
mentum component induced by the tensor interaction. 
On the other hand, when the extra neutrons occupy the 
Pi/2 orbit, these neutrons disturb the tensor correlation 
in "^Hc in some amount, because of the Pauli-blocking 
effect between the extra pi/2 orbit neutron and the Pi/2 
component largely excited from the s-orbit in ^He by the 
tensor interaction. The extra neutrons cannot share the 
spatially compact pi/2 orbit, which is used in ^Hc to pro- 
duce the 2p2h excitations. As a result, the total system 
cannot produce a large enhancement of the tensor en- 
ergy. Therefore, the contributions of tensor interaction 
in each level of He isotopes depends on the neutron con- 
figurations. This state dependence works to produce the 
observed LS splitting energy. The amount of the high 
momentum component from the tensor interaction also 
depends on the neutron configurations and explains the 
kinetic energies of each level. 

IV. SUMMARY 

We have developed a method to describe nuclei with 
bare NN interaction on the basis of the tensor opti- 
mized shell model with the unitary correlation operator 
method, TOSM-I-UCOM. We have treated the tensor in- 
teraction in terms of TOSM, in which the 2p2h states are 
fully optimized to describe the deuteron-like tensor corre- 



lation. The short-range repulsion in the NN interaction 
is treated in the central correlation part of UCOM. We 
have shown the rehability of TOSM-I-UCOM using the 
AV8' interaction to investigate the structures of He iso- 
topes. It is found that the excitation energy spectra are 
nicely reproduced. 

It is found that ^He contains a relatively large amount 
of the pn pair in the pi/2 orbit. The large mixing of pi/2 
orbit in ^He contributes to the production of the splitting 
energy between the LS partners in neutron-rich He iso- 
topes. In ^He, it is found that the 3/2^ state gains more 
of the tensor energy than the l/2~ case. This is due to 
the Pauli-blocking effect between the tensor correlation 
in "^Hc and the motion of the last neutron occupying the 
Pi/2 orbit. In the 3/2^ state of ^He, the enhancement 
of the kinetic energy is observed because of the high mo- 
mentum component brought by the tensor interaction. 
As a result, tensor interaction dynamically produces the 
state dependence in ^He and contributes to the creation 
of the splitting energy between the l/2~ and 3/2~ states 
in ^Hc. The similar mechanism is found in ^~®He be- 
tween the configurations seen in the ground and the ex- 
cited states. The configurations involving manly the P3/2 
orbit of extra neutrons gains the tensor energies and pro- 
duces the large kinetic energies, which are located near 
the ground state region in He isotopes. On the other 
hand, the states involving large Pi/2 components of ex- 
tra neutrons are mostly located in the excited states of 
He isotopes. In those excited states, the mixing of the 
Pi/2 orbit neutrons reduces the enhancement of the ten- 
sor energies. We note that this Pauli-blocking effect of 
extra neutrons with those states in 2p2h excitations due 
to the strong tensor interaction does not exist for the case 
of the Minnesota interaction. 

The enhancement of the kinetic energy due to the pres- 
ence of the tensor interaction indicates the increase of the 
high momentum component in the wave function. Obser- 
vation of the high momentum component experimentally 
in finite nuclei in order to confirm the existence of the 
strong tensor correlation is desired Q . Based on the 
results obtained in TOSM-I-UCOM for He isotopes, it is 
interestin g to include the continuum component of extra 
neutrons [30j and the genuine three-body interaction. 
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